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ABSTRACT 


In ordinary gas dynamics we use assumptions which also agree with 
kinetic theory of gases for small mean free paths of the air molecules. 
The air forces thus calculated have to "be re-examined, if the mean free 
path is comparable with the dimensions of the moving body or even with 
its boundary layer. This case is very difficult to calculate. The con- 
ditions, however, are more simple, if the mean free path is large com- 
pared to the body length, so that the collisions of molecules with each 
other can be neglected compared to the collisions with the body surface. 

In order to study the influence of the large mean free path, calcu- 
lations are first carried out for the case of extreme rarefaction. 
Furthermore, the calculations on this "completely ideal" gas will be 
carried out under consideration of •Maxwell , s velocity distribution and 
under the assumption of certain experimentally established reflection 
laws for the translational and nontranslational molecular degrees of 
freedom. 

The results thus obtained allow us to find, besides the air pressure 
forces perpendicular to the surface, also the friction stresses parallel 
to the surface. 

Their general results are calculated out for two practically 
important cases: for the thin smooth plate and for a projectile— shaped 

body moving axially. 

The mathematical part of the investigation was carried out primarily 
by Dr. Irene Bredt. 


*"Gaskinetik sehr hoher Fluggeschwindigkeiten. " Forschungsbericht 
Nr. 972, May 31, 1938. 
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I. INTRODUCTION 


With the motion of todies at very great atmospheric heights, the ■ 
air can no longer he considered a continuous medium, in the sense of 
flow theory. 

At over 50 kilometers altitude, the mean free path of the. air 
molecules will he of the magnitude of boundary layer thickness and, 
at over 100 kilometers altitude, of the magnitude of the moving tody 
itself . 

The mean free path at greater heights will he definitely greater 
than the tody dimensions of the moving tody, and the especially simple 
conditions of very rarefied "completely ideal" gases are valid where 
the effect of the collisions of the molecules with each other disappears 
compared to the effect of the collisions with the moving tody. 

The air molecules collide then against the moving tody as individual 
particles, independent of each, other, and are reflected with a mechanism 
which deviates more or less from the known Newtonian principle of air. 
drag, as shown ty the results of the kinetics of very rarefied gases. 

The velocity of the tody will te denoted by v (meters per second) 
in the following discussion. 

If we imagine as usual that for the consideration of_.the flew 
process the tody stands still and the medium moves, then v equals the 
uniform undisturbed flow velocity of all the air molecules. 
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The air molecules have also their ideal random thermal motion. 


The individual molecular thermal velocities are distributed com- 
pletely at random in al 3 directions and are of completely arbitrary 
magnitudes, where, the various absolute velocity magnitudes group around 
a most probable value, c (meters per second), according to Maxwell* s 
distribution, which has the ratio y[2: JT to the most applied gas 
kinetic value of "average molecular speed c" (square root of the average 
of the squares of a~l 1 the velocities present). 


According to the known Maxwell speed distribution law, letting 
p(kg s^/fflh) equal the total molecular masB per unit volume, the 
mass dp of those molecules having velocities between c x 
and c x + dc x is: 


p ' 

dp k c-r 

~ = -==■ e 
P \fn c 3 



( 1 ) 


If we consider only this mass dp of an otherwise motion- 
less (v = 0) gas, the molecules of which are moving with the particular 
speed c x in random directions, then the quantity dp of molecules 
striking per second on a unit surface of a motionless plane wall can 
be calculated, imaging that all velocities c x are plotted from the 

center of a sphere with radius c x . The conditions of figure 1 are the 

result of an inclination angle 0 between the wall, normal and the 
velocity direction under consideration, and of the molecular quan- 
tity dp — cT~ which passes through the striped area dJF of the sphere 

^c x 2 * 

surface. Then: 


.0 = 


Prt/2 


dp 


he. 


C T cos 


0=0 


0 <3JT = if C* 
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If, with, the aid .of Maxwell* s equation, we include in our "calcu- 
lation ad 1. molecules with the various possible speeds c x between - 0 
and °°, the total mass p of the molecules colliding per second will 
be: 




P = 


x 

4 


c x=° 


00 



°x=° 




PC 

2\[H 


( 2 ) 


as one can find in any textbook of gas kinetics. The pressure of the 
motionless gas against the stationary wall can be calculated similarly. 

The impulse dip .perpendicular to the wall with which molecules of 
particular speed c^- are striking the wall at an angle f will be: 


din = 


dp 

4°/* 


'rt/2 


Cyf COS 2 0 dF = C x ‘ 


0=0 


and the impulse of all molecules striking the wall: 


P 00 


s- 


c x 2 
— dp 

6 P 


Jc-w — 0 


a _p_ 

3 yit 


c. 

c2 


.4 “ 


''X 

.2 


= £ c 2 (3) 


dc v — t— c 
x 4 


°x=0 


If we double this impulse value because 
always assumed for motionless gases, then we 


of the elastic rebound 
obtain the gas rest pressure 



The value of the total impulse is also interesting, i.e., the sum 
of all single molecular impulses which strike per second on the unit - 
area. 
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The total impulse of the molecular mass dp correspond i ng to a 
particular c x will he: 


di = dpc x = 



It is thus 1.5 times greater than the effective impulse dip 

against the plate. The total impulse i is also greater in the same 
ratio. 


II. AIR FORCES ON THE FRONT SIDE OF A FLAT PLATE 
OBLIQUE TO THE AIR STREAM 

\ 


If we consider again the mass dp of molecules with speeds 
between c x and c x + dc x (almost equal) and if we examine its action 
on a flat plate in an air stream with an angle of attack a, then this 
process can he illustrated by figure 2, if we further assume 
that v sin a < c x . 

The uniform velocity v of the individual molecule combines with 
the ideal random velocity (which can have any space direction) to give 
a resultant , the components of which are: 

perpendicular to the plate: 
parallel to the plate and to v cos a: 
parallel to the plate and perpendicular 
to v cos a: 

and for which the absolute value is therefore: 

w = |/v 2 + 2vc x (sin a cos 0 + cos a sin 0~ cos t|t + c x 2 ) 

From the sphere of all possible directions of q x , a spherical sector 
with the half opening angle cob af = v sin a/c x is excluded, in which 
the speed component v sin a + c x cos 0 is directed away from the 
plate, i.e., the molecules of this 0 range do not strike the 'plate. 
The integration over the velocity directions of all colliding molecules 
is not from 0=0 to & as with the motionless gas, but from 0=0 

to Jt = 


v sin a + c x cos 0 
v cos a + c x sin 0 

c x sin 0 sin T)r 


cos t 
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The .molecular mass, colliding per second against the unit plate 
area with the selected speed c x is therefore: 


dp = 


dp 


ik3 x 2 Jt 


‘n-Hb 


(t sin a, + c x cos t )2c x 2 it sin 0 d0 


JZUO 




c x + 2 t sin a + 


T^sin^a 




For v sin a > c x the integration extends over the whole sphere 
from 0=0 to it and the molecular mass colliding against the plate 
with selected- speed c x is 


dp = 


dp 


hc x 2 it 


pit 


0=0 


(v sin a + c x cos 0)2c x 2 it sin 0 d0 = dpv sin a 
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Both equations naturally give the same value for v sin a = c x . 

With the aid of Maxwell's distribution equation the total molecular 
mass colliding on the unit area per unit time within the speed 
range c x = 0 to » will be: 


n co 


P = 


Cj—v sin a 


K c 


_ . Y 2 sin^aA , 

c x + 2v sin a + ]dp + 

c x 


Hv sin a 


v sin a dp 


c x =° 


p 00 


= - 2 — 
\f*~ 


4p 


c x =v sin a 


(: 


e v 3 2vc. 


sin a + 


v^c 


x . —X . ■ -X eln 2 a ) e 02 d 




Pv sin a 2 




vc. 


sin a e c dc 


Ox-0 




■^si^a 


’ v sin a, c x ^ 


P cl C 2 

___ — \e 

yir ^ 


r— v sin a , 0 v Bin a 
+ \J rt + 2 


e c dc. 


c x =0 


Thus the number of the colliding molecules is known, and for cal- 
culation of the forces acting on the plate, we must now determine what 
impulse the molecular mass under consideration produces in the direc- 
tions in question. 

The impulse perpendicular to the plate, ip and the impulse 
parallel to it i T will be ex amin ed separately. 
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We find for the Impulse per second perpendicular to the plate 

■'it-se 


1. If t sin a < c x : dip = 


dp 


4c x ^it 


/0= 0 


(v sin a + c x cos 0)^dF 


• f * £ B1Xi a ) : 


2. If v sin a > c T : di. 


dp 




(v ein a + c x cos 0)^dF 


x * J0=O 


= —■ cj^^2 + 6 sin^i 


6 * V Cx 2 


“) 


This summation of the impulse components over all possible direc- 
tions yields the total impulse., perpendicular to the plate ^ of the air 
molecules striking the unit area per second with. a speed c x deter- 
mined by dp. 
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The further summation of the impulses over »i T speeds c z -with the 
aid of Maxwell’s distribution equation results in the total impulse i 
acting per second perpendicular to the plate: 




If we set a = 0 then equation ( 5 ) gives the impulse of the motion- 
less gas against the motionless wall, equation (3). 

Also a = 0 gives the impulse of the motionless gas, which will 
not he influenced through uniform motion of the gas mass parallel to 
the flat plate. We find the impulse parallel to the plate in the direc- 
tion of v cos a in the following manner: 

The molecular beam with a particular speed c x and with a particu- 
lar direction (the latter determined by the inclination angle ) between 
the velocity c x and the perpendicular to the plate, and by the angle \|r 
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between the projection c x cos 0 and the direction v cos a) gives the 
molecular mass colliding per second according to figure 2. 


dp = dp 


c x ^sin 0 d0 dijr 


( T 


sin a + c x cos 


0 ) 


Hie velocity component of this beam parallel to the plate is: 
v cos a + c x sin 0 cos i(r 

and the impulse of the beam with a particular c XJ 0, and will be 

dp s - — 0 - sin a + c x cos 0^v cos a + c^- sin 0 cos 

If one integrates over all i|r and f, one obtains the impulse of 
the beam with a particular c x . 

For v sin a < c x : 


P2n 


di 


dp 




Prt-^ 


Ji|r=0 


0=0 


^v sin a+ c x cos 0) ^v cos a+c x sin 0 cos ^sin 0 d0 diK 


dp „ 2 
6 


( 3 v3 p v 2 3 v \ 

I — sin^a cos a + 3 sin a cos a + — — cos al 

\2 Cx 3 c x 2 2 c x . ) 


and for v sin a > c x : 




2jt 




0=0 


^v sin a + c x cos 0)(v cos a+c x sin 0 cos lf/^sin 0 d0 d-f 


= dp v 2 sin a cos a 
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The integration over al 1 c x with the aid of Maxwell* s distribution 
equation gives fina l ly the total impulse in the required direction 
parallel to the plate: 


Pv sin a 


i T = 


dpv 2 sin a cos a + 


c x=° 




c-^v sin a 


— c x 2 (^ — =■ sin 2 a cos a 
6 V Cx - 3 


+ 3 



sin 


a cos a + 


1 -Z_ 

2 c x 


cos a 




v 2 sin 2 a 


= P 


2 \f* 


vc cos a e 


/ 

"*v sin a — ^ \ 
c 2 \ 

e dc 

( — v^sin a cos ajll -i 2 — 

\ 2 A c ^ 

/_ 


( 6 ) 


We can start out from, the impulse of the gas stream, given "by 
equations ( 5 ) and (6), (perpendicular and parallel to the plate) in 
order to calculate the forces produced by the air on the front side of 
the plate oblique to the air streamy including the force perpendicular 
to the plate (normal pressure p) and the force parallel (friction 
stress t ) . 


For this calculation, we have to make same assumptions an the trans- 
fer of this impulse to the plate and on the change of the kinetic trans- 
lational energy of the moJLecules into other energy forms for which pre- 
sent gas kinetics furnish only a partial basis. 

If we first assume monatomic gases, so that inner degrees of free- 
dom for energy absorbance do not exist, and further assume that the 
struck molecules of the wall are in such a temperature condition' that 
they also cannot take over any energy from the colliding molecules. 
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then the molecules have to leave the wall again with the same speed with 
which they arrived. The collision is thus completely elastic and we have 
only to derive the direction of reflection. 

Gas kinetics distinguishes two different possibilities for this: 

Mirror Reflection where the assumption is made (following Newton) 
that the angle of incidence and the angle of reflection are equal and 
both beams are in the same perpendicular plane. 

Diffuse Reflection where it is assumed (following Khudseni) that 
the reflection direction is not at all dependent on the direction of 
the impinging beam and is completely diffuse, i.e., that the colliding 
molecules first submerge in the wall surfaces, then after a finite time 
of "adherence" leave again, in a completely arbitrary direction. 

This last hypothesis is generally accepted in flow theory, where 
the adherence of the frictional boundary layer on the surface is 
explained by diffuse recoil of the molecules. 

Xh the case of a motionless gas (v = 0) both assumptions lead to 
the same distribution of rebound molecules and thus to the same forces 
bn the wall, as the striking molecules are in completely random direc- 
tions and this then is also true for the rebound molecules under both 
assumptions. ~ 

In the case of a gas in motion, the two assumptions lead to very 
different air forces. 

With elastic "mirror" reflections, the impulse 1 T of the gas 
flow parallel to the wall stays unchanged. Shear stresses on the front 
side of the plate t t are not transferred to the wall. Okie friction 
forces are zero. 


T 

V 


0 


(7) 


The impulse i^ of the impinging molecules normal to the wall is 
destroyed completely and an equal but opposite impulse is produced by 


^Knudsen, M. : Annalen der Physik, Yol. 28, p. 75, (1909); Vol. 28, 

p. ii4, (1909); Yol. 28, p. 999 , (1909); Vol, 31. pp. 205, 633, (1910); 
Vol. 35 j p. 389, (1911); Vol. 34, P- 593, (1911); Vol. 48, p. 1113, 
(1915); Vol. 50, p. 472, (1916); and Yol. 83, p. 797, (1927). 
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the completely elastic recoil. The pressure on the wall caused "by this 
process thus equals twice the impulse ip: 


Pt = 2i p 


( 8 ) 


With diffuse-elastic reflection, the impulse i,. of the gas stream 

parallel to the wall will he giTen up entirely to the wall and the fric- 
tion stress equals Ij. : 


T v = V 


( 9 ) 


The wall normal impulse ip of the arriving molecules is destroyed 
again, whereby is created, a partial pressure p.j_ = ip. 

The second, part of the pressure , due to the diffuse— elastic recoil , 
has to he investigated more closely. To imitate the real process, an 
impulse value of the magnitude of the complete impulse i of the beam 
striking per unit plate area is distributed evenly on a hemisphere as 
if all gas particles started from the center of this hemisphere, and 
finally the resultant of this impulse distribution perpendicular to the 
plate is ascertained. 

The total impulse i of the arriving beam is derived, according 
to the proceeding impulse calculations, in the following manner: 

The molecular mass striking per unit time on the known area 
section df = c x 2 sin 0 d0 dty selected for a certain c x , 0> and ijf is 

dp = dp fc x 2 sin 0 d 0 dty Ik c x ^ it) (v sin a + c x cos 0) 


The effective speed of this ray is 
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and tits impulse per second' thus: 


dpw = dp 3in -j~ — (y sin a 


+ c x cos 



+ 2vc x (sin a cos 0 + cos a sin 0 cos 


+ c. 


This impulse integrated over all c x , 0, and ^ will give finally the 
total impulse of the team. The actual carrying out of this integration 
is so difficult that the impulse will evaluate hy successive approxi- 
mation. It is started with the vector sum of the impulse resultants i 

and I T , perpendicular and parallel to the plate. . 


r 2 = i. 


V 


This impulse resultant is smaller than the total impulse. It is 
found in connection with equation ( 3 ) that the total impulse of a gas 
at rest (v = 0) is 1.5 times the impulse resultant. Total impulse and 
resultant are equal for uniformity flowing gas without heat 
motion (c x = 0) . For conditions lying in "between, we assume a constant 

relationship for the factor with which the impulse resultant has to he 
multiplied to obtain the total impulse i. For instance: 2 


i = 


1.5 + g- sin. a + 1.5(iy) 2 - 

1 + J em a + 1 


'This total impulse, according to our assumption, is now .considered as 
the completely uniform impulse radiation per area unit in ai i directions 
. outward from the surface. 


2 

Interpolation by Prof. Busemann, Braunschweig. 
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The "beam, pressure pg perpendicular to the surface is then equal, 
to i/ 2 , as shown "by a simple integration over all normal components. 

Thus, the pressure vertical to the wall in the case of diffuse- 
elastic reflection is: 


Pv = 1 + I/? 


( 10 ) 


In order to get a first view on the numerous, conditions of the air 
forces Just found, figure 3 shows the relationship "between the pres- 
sure p or the shear force t and the dynamic pressure q = £■ v 2 for 

a hypothetical atmosphere of monatomic hydrogen with t = 0° C tempera- 
ture (c sb 212k m/s) and for flight speeds up to v = 8000 m/s, for 
either mirrorlike or diffuse recoil. 

It Is seen how different the air forces can be according to the 
assumptions made: mirrorlike or diffuse. 

In gas kinetics an attempt is made to approach the real conditions 
by assumi n g that the reflection for a fraction f of all striking 
molecules is diffuse, while the remainder (l — f ) will be repelled 
mirrorlike. The fraction of diffuse reflections depends on the kind 
of striking molecules and particularly on the material, surface conditions, 
and temperature of the struck wall. 

According to numerous measurements3, the plate can be considered as 
completely rough under conditions usually prevailing in flight technique, 
i.e., .the mirror reflected part (l — f) Is negligibly small, so the 
reflection will be almost completely diffuse. 

An experimentally obtained dependence of f on the angle of attack, 
such that the reflection will be more mirrorlike with smaller angle of 
attack, is according to previous measurements of flight relations, too 
inslgnif leant to be considered. 

Knauer and SternA assume from optical analogies that the angle of 
attack at which mirror reflection begins is such that the surface roughness 

3for example, Karl Jel l i n ek, Lehrbuch der physikalischen Chemle, 

Yol. 1, p. 270, 1928. 

^Knauer, F.., and Stem, 0.: Zs. f. Phys. Yol. 50, pp. 766 , 799 

(1929). 
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height, projected an the beam, must be smaller than De Broglie's — 

wavelength X of the molecular beam. 

With a wave length of 10~ ® cm and a roughness height of 1 0“5 to 10~^, 

one obtains for the angle of attach which is of interest sin a = X/h = 10~3 
i.e., an angle range of a few minutes, which is insignificant in flight 
technique. 

With regard to the reflection direction, we assume in the following 
analysis that f = 1, i.e., completely diffuse recoil. 

So far, for the reflection speed, perfect elasticity of the recoil 
was assumed, which means individual recoil speed is equal to the 
colliding speed. 

The struck wall will, in fact, be much colder than the gas molecule 
temperature after its submergence in the wall surface (and so after its 
complete braking an the plate to the velocity v), so that we have to 
assume heat transmission to the wall molecules from the colliding 
molecule which remains a finite time in the plate surface. 

Figures k and 5 show the internal energies U for molecular 
nitrogen or hydrogen and their combination from the individual degrees 
of freedom of the molecular motions as a function of the temperature, 
starting from an inte r nal energy U 0 of the gas at rest corresponding 
to a temperature of 0°C, with the other internal energy values equal 
to the kinetic energy corresponding to v. 

For this, the relation between U and v is stated 
as U = U 0 + Av^ / 2g . The graph goes up to a U = 8000 kcal/kg, corres- 
ponding to a flight speed range up to v = 8000 m/s. The specific heat 
at constant volume c T was calculated under the usual assumptions on 

energy absorption by translation, rotation, and oscillation of the 
molecules (the latter according to Planck's formula) after the collision. 

We see from both figures that, especially for the very high 

temperatures correspond to the high flight speeds. A complete tempera- 
ture equalization of the colliding molecule to the wall temperature 
would he equal to a total annihilation of the recoil speed, or an 
almost completely inelastic collision. It should be observed that the 
wall accommodates itself in a Bhort time to the temperature of the 
colliding molecules, because of the very small heat capacity of the 
thin metal wadis of the moving body. 

The molecule mass, colliding on the oblique unit surface per second 
at very high flight speeds is, for exanq>le, about equal to pv sin a. 
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and tins the arriving energy S = TJppv sin a. With pp = 7 = 10 “ ^ kg/m. 3 , 
a = 7 °, and v = 8000 m/s the value of the energy brought in is 

1 = 7.8 kcal/m 2 sec, almost independent of the composition of the atmos- 
phere. If the accommodation coefficient of the arriving gas molecules 
is one, then the wall would obtain this energy in the form of heat and 
this heat quantity should "be given away 'by radiation, where a temperature 
increase AT of approximately 580° is necessary for black body radia- 
tion, i.e., the plate stays in fact pretty cold compared to the colliding 
molecules, and a very intensive, lasting energy delivery by the colliding 
molecules is. out of the question. According to existing measurements 5 , 
this temperature equalization is not 100 percent, however, an accommodation 
coefficient of 30 percent was found under certain conditions,' i.e., the 
reflected gas mass contains still 70 percent of its internal energy IT 
which it possessed at the moment of collision. 

The reflection velocity for a monatomic atmosphere is established 
this way. 

The remaining internal energy of a molecular atmosphere will dis- 
tribute itself quite differently over the existing degrees of freedom, of 
the reflected molecule than assumed for the colliding energy, which con— . 
sisted primarily of kinetic energy l /2 Av 2 /g, and only in small, measure 
of the internal energy U 0 of the gas at rest, which latter was distri- 
buted evenly over all degrees of freedom. 

For a diatomic molecule with three translational and two rotational 
degrees of freedom, the individual shares of U 0 , for an "average" 
velocity c, are for each kilogram of gas 3/6 Ac^/g for the three 
translational degrees and 2/6 Ac^/g for the two rotational degrees of 
freedom of the molecules. 

On collision, all degrees of freedom will take part in the energy 
distribution change, and it can be assumed for further estimation of 
the diatomic molecule between the perfectly elastic and the perfectly 
inelastic collision, conditions, for Instance, that the total energy 

A/gQj- v 2 ’+ 3 . C 2 + ir 2 ^ = A/g^- v 2 + ^ ff 2 ^ distributes itself on the. 
average evenly over all these five degrees of freedom. 

It can be taken from figures 4 and 5 that very high temperatures 
are associated with the high colliding molecular speeds at which another 
motions,] degree of freedom is excited, that of mutual molecular oscillation 

''Vien-Harms, Handbuch der Experimentalphysik, Tol. VIIl/2, p. 638, 

1929 . 
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Hie known Boltzmann *s equilization rule cn the energies is not of 
value for these oscillational degrees of freedom. 

While the three translational and the two rotational degrees of - 
freedom of a diatomic molecule have each the same energy admission 

TJ X = - ART = - Ac 2 /g (kcal/kg) 

2 6 


that Is together 


U (trans. + rot,) 2 ^ = 6 


the energy admission of the oscillational degree of freedom TJ 0 at low 

temperatures Is practically zero and increases at higher temperatures 
according to Planck's equation: 


TJ 


8 


ASQ 2 A -2 £ 

e? /T _ 1 " 6 « t („«A _ i) 


approaching the limiting value, valid for high temperatures, of full 
exitation of the oscillational degree of freedom. 


TJ„ = - ABT = 


6 g — . 


In the last equation: 

A = the mechanical equivalent of heat, 1/^27 kcal/kg 
E = the individual gas constant m/° 

T = the absolute temperature °K 

$ = a temperature characteristic for each material, which is, for 

example, for nitrogen Eg = 3350 °£> for hydrogen Hg = 6l00°K. 

The temperatures of the colliding molecules at high flight speeds 
are so high (according to figs. 4 and 5) that the molecular gas here 
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already dissociates strongly into itB atoms under normal equilibrium 
conditions . 

The transformation of the gases hydrogen and nitrogen, which are of 
importance in the higher atmospheric layers, into their monatomic, active 
modification belongs to the most energetic endothermic chemical processes 
which are known (Hg = 2H — 51300 kcal/kgj Ng = 2W — 6050 kcal/kg), and 

the dissociation (if it actually occurs) would absorb extraordinary amounts 
of energy and would make the collision almost completely inelastic. 

So far it has not yet been proven by experiments that these 
molecules really dissociate an collision against a fixed wall at the 
speeds here considered. 

However the results of known tests with electrons colliding 
against the molecules of very rarefied Hg or Ng gases let us 

guess that the energy of the collision with a molecular speed up 
to v = 8000 m/s is not sufficient to disturb the molecular bond. 

With the electrons colliding against Wg or Hg molecules dis- 
sociations are observed^ only when the energy of the colliding electron 
was several times the dissociation energy of the struck molecule. 

This transferred on our case would yield colliding speeds of 
over v = 10,000 m/s for a nitrogen atmosphere or over v = 35,000 m/s 
for a hydrogen atmosphere, which lies outside of the range of our 
investigations . 

For the calculation of the forces on a plate oblique to the air- 
stream, we shall therefore not assume the dissociation of the colliding 
molecules . 

The degree of elasticity of the recoil will only be derived from 
the energy distribution of the wall molecules and of the proper rotational 
and oscillational degrees of freedom of the colliding molecule. 

This degree of elasticity, i.e., the ratio of the molecular 
reflection velocity, when energy division occurs, to the reflection 
velocity when ho energy division occurs, is estimated as follows: 


Wien-Harms ; Handbuch der Experimentalphysik, Tol. YIU/l, 

pp. 704, 706 (1929). 
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First, corresponding to tlie measured accomodation, coefficient, 
let 30 percent of the collision energy corresponding to v "be ccraamnicated 
to the v ail molecules. The remaining collision energy, 

A/ y 2 + 1^2 + |c 2 J= A/g^il Y 2 '+ | 5 2 ) 

should he distributed over the three translational degrees, the two 
rotati onal degrees and the oscillations! degree of freedom, evenly and 
according to the exitation degree so that each translational, degree of 
freedom has the following energy: 


The condition of motion of the diffusely reflected molecules is the same 
as assumed with equation (3), only there the energy content of a trans- 
lational degree of freedom was A/g -g C 2 . Here the internal energy 

belonging to the three translational degrees of freedom is given 
as U = 3E r and the translational speed is now: 


5 + 22 

*(e‘ /T - l) 




instead of: c r + c 2 " when no energy distribution takes place 

between the vail and the translational degrees of freedom. (The latter 
value comes from the energy balance: 
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The recoil impulses should fall off approximately like speeds 
obtained from energetic considerations , so that: 



(H) 



This degree of elasticity is drawn in figure 6 for nitrogen and 
hydrogen as a function of the flight speed t, using figure h or 5 for 
the relation between v and T. 

The pressure vertical to the wall on the plate oblique to the 
air stream, in the case of diffuse— inelastic reflection, can be obtained 
from equation 11: 


Pv = 



Corresponding to figure 3.» figure 7 shows the relation between the 
pressure p or the shear r and the dynamic pressure q for hydrogen 
(o *3 1508 m/s) and nitrogen (c = ho6 m/s) for air angle of attack a = h° 
• and for flight speeds to v = 8000 m/s. 

In figure 7 the two impulse contributions due to collision (ip) 
and recoil (ei/2) are separated. 
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Hie collision impulse could be given, without objection from purely 
mechanical relationships. 

t 

The reflection impulse is estimated on the basis of a series of 
rather arbitrary assumptions on direction and speed of the reflection. 

Further inf ormation, proceeding out of .the conjectures presented 
here, cn the retention of single air molecules after collision with 
the w all with high, speed should be found by experiment. Such tests may 
be joined with the well-known, molecular beam investigations, where the 
usual beam speeds are to be increased greatly by correspondingly greater 
energy delivery to the molecules under investigation. One can thus 
obtain a type ‘of wind tunnel investigation in which single molecules 
fly with extremely high velocities, and the effects of their collision 
with a solid body can be observed. 

With the help of very rapid molecular beams, a number of questions 
should be cleared up experimentally: how often a reflection of the 

molecules from the struck wall actually occurs, what factors change the 
completely elastic collision into a more or less Inelastic collision 
through the transformation of translational energy into other forms of 
energy (i.e., rotational, oscillational or dissociaticnal energy of 
the gas molecules or the wall molecules), and what direction law the 
fi n al reflection follows, whether mirrorlike reflection, or predominantly 
diffuse reflection, or following a different law. 

These investigations can use to advantage De Broglie* s analogy 
between molecular beams and x— rays. 


IH. AIR FORCES ON TEE FRONT SIDE OF A FLAT PLATE 
VERTICAL TO THE AIR STREAM 


Vertical flow against the flat plate represents the limiting 

of oblique flow. 

The relations obtained in section 2 are preferably discussed her® 
for this special case. 
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The total molecule mass, colliding per unit area per unit time 1b 
given "by equation 4: 



The complete integration of equation 4(a) (in contrast to equation 4) 
is possible by series development, because one can put v/c»l for the 
high flight speeds under consideration, -which -was not possible for 
v sin a/c in equation 4. 
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Hie total Impulse of this gas mass striking vertically against the 
plate per unit time and unit area is given "by equation 5: 




pc‘ 







1 £ 1x3 £ 

2 -y 2 2 2 v^ - 



(5a) 


If , as assumed at the start of the integration v/o » 1, the terms 

_ S_ 

of the equation multiplied "by e c 2 can he neglected compared with 
the first two terms of* the equation. 


For instance, the influence of the variety of absolute molecular 
velocities (following Maxwell's distribution), which Is presented by 
these higher terms, is less than 0.3 percent. of the value of the first 
two terms, when v/c is 2 or more. 

Therefore, for v/c > 2, the impulse can be calculated as if all 
the molecules had the same velocity cj thus, considering only the 
first two terms of equation 5(a), 


ip = pv 2 + | pc 2 (5b) 


■Under the assumption of a completely inelastic collision, this 
impulse Is equal to the required pressure on the plate: 


P = 


ip = pv 2 + | pc 2 


(10a) 
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If It is assumed that the air molecules have not lost their Telocity 
after the collision -with the plate, hut that they rehound perfectly 
elastically and "mirror— like," then the impulse given to the plate is 
simply doubled, and equation 10(a) for elastic mirrorlike collision 
is -written: 


p = 2ip = 2^pv2 + | c 2 ) 


(lOh) 


The influence of the random, molecular motion, represented by the 

second term, is at ^ = 2 approximately 12.5 percent of the pure 

c 

Newtonian air force, which is represented by the first term. However, 

the influence decreases, from ^ = 5 on, to under 2 percent of the 

c 

Newtonian value, contrary to the oblique conditions with small angles 
of attack, where the influence of the molecular velocity is still great 
even with high flight speeds. 

If it is assumed that the molecule reflection is completely elastic 
but diffuse, then the pressure on the plate decreases to a value: 



+ 0.75 


2 v 

1 + ^ - sin a + 


1 + - sin 
o 





= plv 2 + |c 2 ) 



0.75 


1+ |Isin a+ (l) 2 


1 + — sin a + —(—') J ’ 

o 2 \ o / / 


(10c) 
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If finally the molecule reflection "becomes not only diffuse "but also _ 
partly inelastic in the sense of equation 12, then the pressure on the 
plate reduces again to; 



The graphs in figures 8 and 9 show the pressure conditions according 
to equations 10(a) to 10(d) versus the dynamic pressure of nitrogen or 
hydrogen. 

For the calculation of the pressure of a motionless gas on the walls 
of the gas container which has the same temperature as the gas kinetics 
assumes completely elastic collisions for the normal range of molecular 
speeds, i.e., the molecule keeps the same translational energy, an an 
average, as It had "before the shock. It is not important to know if 
the reflection is mirrorlike or diffuse, because both assumptions lead to 
the same molecular picture for the gas at rest. 

If the mplecules are polyatomic then the distribution of the total 
energy to each possible degree of freedom is already uniform before the 
shock j and this uniform distribution need not change after the shock. 

In equation 10(a) (for completely inelastic collisions) for the air 
pressure against a plate vertical to the air stream p = pv^ + £■ c^ 

the first term pv^ corresponds to the dynamic pressure of molecules having 
no random motion against the plate (Newton), while the second term ^ pc^ 

corresponds exactly to the pressure of the static atmospheric air. 
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However, this explanation of the individual terms Is only for mally 
correct , since the resting air pressure is calculated assuming elastic 
molecular collisions. With the inelastic collision, the decrease in 
the "stopping" pressure due to loss of the recoil impulse will he off- 
set by the mixed term, in the square of the sum of the two speeds (y and c), 
in the velocity range under consideration. 


This condition can he recognized more clearly from equation 10(h) 
for the completely elastic collision p = 2pv^ + pc^ where after sub— 

1 9 

traction of the static air pressure — pc c there remains a pressure 

P 1 o p 

of 2pv + pc% which contains besides the Newtonian term, 2pv , 

^ 12 

also an additional term of — pc , which reflects the effeot of the 
mixed term. 


No resulting impulse is given parallel to the vertical plate 
because of the symmetry of the total system, i . e . , friction forces are 
transferred in the plate plane, hut the sum of these forces outward 
is zero. 


IV. AIR FORCES ON THE BACK SIDE OF A FIAT PIATE 


If the mean free path of a molecule iB small compared to the 
dimensions of an empty space into which the gas is flowing, then the 
flow— in speed can he greater than the most probable molecule speed. In 
the flow of diatomic gases into a complete vacuum, the directed 
velocity, of the total flowing mass can surpass the probable 

molecule speed c by a factor of about 1 . 87 , according to the laws of 
gas dynamics. 

If, however, the molecular mean free path is comparable to the 
empty spare dimensions or even greater than these as assumed here, then 
the number of molecule collisions behind the rapidly moving plate during 
the flow— in is not sufficient to produce the mentioned acceleration, 
and the molecules move with their usual speed c into the empty space 
behind the plate. 

According to figure 10, collisions between the air molecules and 
the back side of the plate cannot take place, i.e., the pressure against 
the "suction side" of the plate must have become zero, as soon as 

> ° 
v = 

sin a 
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This "border line is, however, strongly blurred because of Maxwell’s 
distribution. 

The effective forces against the back side of the plate can be 
derived by the same process which led to equations 4 to 12, where now . 

the uniform velocity v is directed away from, the plate at the angle a, 

while before it was directed toward the plate (fig. II). 

It is assumed first that v sin a <c 2 * 

The uniform .velocity v of the individual molecule combines with 
the ideal random velocity c x of the molecule (which can have any space 

direction) to a resultant whose components are: 

perpendicular to the plate: c x cos 0 — v sin a 

parallel to the plate and to v cos a: v cos a + c^. sin 0 sin if 

parallel to the plate and perpendicular to . v cos cm c x sin 0 sin i|r 

From the sphere of all possible directions of c^, a spherical sector 
with the half angle cos 3C = v sin — is taken, inside of which the ' 

°x 

velocity component v sin a — c x cos 0 is directed towards the plate, so 
that the molecules of this 0-range do in fact collide against the plate. 
For all 0 > the resulting molecular velocity is directed away from the 
plate. Therefore no collision with the plate takes place. 

The integration over all colliding directions is extended from 0=0 
to 0 

The molecular mass with the chosen speed o x colliding per second, 
against the unit area of the plate Is therefore: 


dp = 


dp 


4c/ rt 


pae 


U0=o 


(c z cos 0 - v sin a )dF = v-(px ~ 2v sin a + v^sl^a/cx) 


For v sin a > 0 ^- there are no collisions with the plate, so that 
this case need not be treated. 
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The total molecule mass colliding on the unit area per unit time 
is obtained with the aid of Maxwell * s distribution equation for the 
velocity range c s = v sin a to » : 


P = 


Jcx=v sin a 


^c x — 2v sin a + v2sin2a/c x ^dp 


P 


fc T =v sin a 


x 

v^c x \ c e 

2 — 9- sin a + sin^ale dc„ 

3 c 3 c 3 / x 


„ TCx 


v^sin^a 


0“ 


P o 
— — — e 
yfx 2 


— x- v sin a/c 

>P 


°x 

e dc. 


Jc^v sin a 


_ o 
2 


v^sin^a 


+ 2v sin a/c 


Pv sin a _ 

e dc. 


c z =0 


2 


(13) 
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Similarly for the impulse vertical to the plate: 


di-n = — I . ( Cx cob $ - v sin <x) 2 dF 
*'*■!** 




■ f °4 - 


3 ~ sin a +• 3 ■— sin 2 a - ~ sin 3 , 






*P" 


lox= V sin a 


6 ° x2 ( X “ 


3 c sin - a + 3 sin 2 a - 3in3a)dp 
z Cx 3 1 


P 


v sin a 


- i v^in a + i cv sin a - - c 2 ){ 

i 3 3 2 / 


■v^sin^ 


+ (| ° + ■ v2sin2a )|\^ 


Pv sin a 5- 

c d 

® dCj;^ 

c=0 


(lfc) 
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Finally, for tie impulse parallel to the plate: 


ox 


clir = 


_ d-P 


lk5TT 2 rt 


fo-y cos 9 — v sin <xjy cos a dF 


J(^=0 


dp/ 

xk 


v cos a — ^v^sin a cos a + — sin^a 

°x 


cos 


a) 


poo 


i T - 


v cos a — 2v 2 sin a cos a + ~~ sin 2 a cos a,^ dp 


(J3j=y sin a 


P 1 c 

— T cos % e 


Y^sit^a 


y sin a 

- c' 


uCj^Y sin a 


dc. 


Xfif 


y cos a 


Y^sin^a 

— rs— 


t* 


y sin a/ 


Pv Bin a 2 


Cx=0 


_ m 

2“ 2 


e o~ dc 


(15) 


The impulse resultant vertical to the wall as a consequence of the 
elastic-diffuse molecular recoil is from the total impulse i of the 
colliding molecules: 


i 


1.5 


1 + | | Bto a + ( 

1 + f .in a + §( if 



+ 


It 2 
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Equation (10) again holds for the total _pres sure against the 
suction side of the plate with elastic-diffuse recoil. (Translator’s 
note: Formula missing in original German report.) While for the 

total shear stress on the suction side of the plate, equation 9 is used: 

T r = V - — 

The influence of a certain inelasticity of recoil can be estimated, 
particularly for small angles of attach by the same procedure which led 
to equation (ll), according to which the degree of inelasticity can also 
be specified. 

For the total pressure against the suction side, equation (12) is 
•valid: 


Pr = ip + 6 | 


Corresponding to figure 7, the graphs in figure 12 show the relation 
between the pressure p or the shear t and the dynamic pressure q 
for hydrogen and nitrogen at an angle of attack a = 4° and flight 
speeds up to v = 8000 m/s. 


Y. APPLICATION EXAMPLES 


With the help of the previously mentioned relations, it is possible 
to estimate all the air forces acting on the surfaces of._a flying body 
of any shape, which is moving at flight altitudes of over 100 km with 
speeds’ between about 2000 m/s and 8000 m/s, if definite assumptions are 
made on the composition of the air at this height. ' ' " ' ‘ 

The air forces were differentiated into those which act perpendicular 
to the surface under observation (pressures), and those which act parallel 
to the surface (friction). 

The pressure stresses as well as the shear stresses were found to 
be a function only of the angle of attack and the flight speed, for a 
particular gas. 

In figures 13 and it is shown this dependence of the air forces on 
all possible angles of attack and on flight- speeds between v =■ 2000 m/s 
and v - 8000 m/s for an atmosphere of molecular hydrogen. _ 
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It is to be noted in figure 13 that the air pressure vertical to 
the plate strongly increases with increasing velocity, even with an 
angle of attack cc, = 0, if the molecular recoil is diffuse. 

This representation can be used as a basis for the calculation of 
air force coefficients for certain flight bodies in hydrogen, treating 
each flat surface section separately, with its own angle of attack, or, 
if the body surface is curved, analyzing it into a great number of 
small areas with individual angles of attack (flat areas or symmetrical 
cone areas), and then investigating these. 

As the simplest example, the flat, infinitely thin plate will be 
treated first. The visual symbols are 

A lift 

W drag 

3P wing surface 

and the air force coefficients are: 



and the glide ratio: ■ 




NACA TM 1270 


34 


In figure 15 are drawn the lift coefficients and in figure 1 6 the 
glide ratios of the flat, thin plate according to the above relations. 


On account of the extraordinarily great shear forces -very had 

Cm . 

glide ratios result , which are approximately s = — =1.9 at 2000 m/s 

c a 

with the most favorable angle of attach and which got -worse at higher 
speeds, for example, at 8000 m/s, e = about 2.7. 


The most favorable angles of attack are compar at ively great at 
small speeds, i.e., at v = 2000 m/s, a = 25° approximately, and decrease 
with increasing speed to about 7° at v = 8000 m/s. 


Similarly to the infinitely thin plate, high speed profiles of 
finite thickness can also be calculated, i.e., wedge-shaped and 
lenticular airfoils. Their air force coefficients hardly deviate from 
those of the smooth plate, if they are of moderate thickness.. 


In general, the wings investigated here in the gas kinetics flow 
range behave worse than in the gas dynaml.cs range, where already the 
glide ratios are worse than in the usual aerodynamical flow region. 

The full effect of this unfortunate behavior will be corrected to 
some extent by a flight technique such that at the high flight velocities 
under consideration, inertial forces are developed by the concave down- 
ward flight path, which support the wing. 

Figure 17 treats the question of how great the air drag is in the 
gas kinetics flow range for a body of rotation (projectile form) moving 
axially, with an ogival nose of three calibers radius and cylindrical 
body, end how far the air drag can be improved by a truncated cone bevel 
at the end of the missile. 

These questions can be easily answered with .the aid of figures 13 
and 14 if the ogive is divided into a large number of truncated cones, 
each of which represents a surface with a definite angle of attack. 

The extraordinary value of the drag coefficient is again striking; 
it can be traced to the very great friction forces in the extremely 
rarefied air. 

A noticeable improvement of the drag coefficient could be obtained 
by beveling the end of the projectile; the improvement is about 7 per- 
cent of the origi n a l value. 

Somewhat more tediously but in basically the same manner, the air 
forces on a projectile, airship, etc., at an oblique angle of attack 
can be determined, using figures 13 and l4. 
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YI. SOMMART 


Tli© air forces on "bodies of arbitrary shape are investigated when 
the bodies move with speeds of 2000 to 3000 m/s in such thin air that 
the mean free path of the air molecules is greater than the dimensions 
of the moving body. 

The air pressure acting perpendicular to the body surface, as well 
as the friction forces acting parallel to the surface, are derived -with 
the aid of the calculation procedure of gas kinetics for surfaces facing 
both toward and away from the air stream at any angle. 

The air forces for an atmosphere of definite composition (molecular 
hydrogen) are calculated as a function of the flight velocity at all 
possible angles of attack of a surface and shown in graphs. 

Thereby the friction stresses between air and body surface prove 
to be of the same magnitude as the dynamic pressure and as the air 
pressures vertical to the body surface, i.e., 300 times greater than in 
the aerodynamic flow range. 

The application of the general calculation results to particular 
technically important oases, like thin airfoils and projectile shapes, 
results in extraordinarily high air drag coefficients and poor glide 
ratios even for the theoretically best wing sections 


Translated by Bureau of Aeronautics 
Technical Information 
Navy Department 
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Figure 1.- Velocity vectors of the thermal motion of molecules of a motionless 
gas and their position relative to a fixed boundary wall. 
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Figure 3.- Air pressures p and shear stresses t on the front side of a flat 
plate at 4° angle of attack in an atmosphere of atomic hydrogen under the 
assumption of elastic diffuse or mirror -like recoil of the atoms from the 
wall. 
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Figure 4.- Colliding speed and associated internal energy of molecular nitrogen 
in relation to the colliding temperature of the gas. 
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Figure 5.- Colliding speed and associated internal energy of molecular 
hydrogen in relation to the colliding temperature of the gas. 
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Figure 8.- Air pressure p on the wall vertical to a stream of molecular 
hydrogen, under various assumptions.on the collision process. 


NACA TM 12T0 


43 



Figure 9.- Air pressure p on the wall vertical to a stream of molecular 
nitrogen, under various assumptions on the collision process. 



Figure 10.- Velocity vectors of the molecular motion for collision on the 

back side of the flat plate. 
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Figure 12.- Air pressures p and shear stress t on the back side of a flat 
plate at 4° angle of attack in an atmosphere of molecular hydrogen or 
nitrogen, under the assumption of diffuse and semielastic molecular recoil 
from the wall. 
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Figure IS.- Coefficient p/q of the air pressure vertical to the plate for 
all angles of attack, and for flight speeds between 2000 m/s arid 
8000 m/s in atmosphere of molecular hydrogen. 
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CYLINDRICAL MISSILE BOTTOM 
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Figure 17.- Coefficients of the pressure drag, friction drag, and total drag 
for a projectile -shaped body of rotation, with different missile bottoms'. 
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